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Abstract 

Flavor singlet and non-singlet axial Ward identities are investigated using the Wilson 
formulation of lattice QCD with Clover 0(a)-improvement, which breaks explicitly 
chiral symmetry. The matching at one-loop order of all the relevant renormalization 
constants with the continuum MS scheme is presented. Our calculations include: 1) 
the contributions arising from the Clover term of the action; 2) the complete mixing 
of the gluon operator GG with the divergence of the singlet axial current; 3) the use 
of both local and extended definitions of the fermionic bilinear operators. A definition 
of the gluon operator GG on the lattice outside the chiral limit is proposed. Our 
definition takes into account the possible power-divergent mixing with the pseudoscalar 
density, generated by the breaking of chiral symmetry. A non-perturbative procedure 
for the evaluation of such mixing constant is worked out. Finally, the renormalization 
properties of the composite insertion of the topological charge operator / d^x GG{x) 
relevant for the lattice calculation of the neutron electric dipole moment, induced by 
the strong CP-violating term of the QCD Lagrangian, are discussed. 
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1 Introduction 



The investigation of the axial anomaly is a powerful tool to unravel the structure of the 
QCD vacuum. Since the latter is highly non-trivial, the conservation of classical currents 
may be spoiled by quantum fluctuations of the vacuum. When this happens an anomaly 
is formed as in case of the fiavor-singlet axial vector current. The axial anomaly can have 
deep consequences on various physical observables, like the large mass of the rj' meson or the 
smallness of the fiavor-singlet nucleon axial coupling constant, which is commonly known as 
the proton spin crisis. 

The connection between the topological structure of the QCD vacuum and the axial 
anomaly has been elucidated by 't Hooft IJ. A very interesting case where such an inter- 
relation shows up, is provided by the neutron electric dipole moment [EDM) generated by 
the so-called ^-term of the QCD Lagrangian [2], which in the Euclidean space is given by 

where is the gluon field strength, c the color octet index (c = 1, ...,8) and 6 a dimen- 
sionless parameter. The 6'-term breaks both parity and time reversal symmetries and 
therefore it can generate a non-vanishing value of the neutron EDM^. Available estimates 
of the relevant matrix element however are based on phenomenological models, as the MIT 
bag model of Ref. |3] or as the effective ttN chiral Lagrangian of Ref. 0. Estimates relying 
on non-perturbative methods based on the fundamental theory, like lattice QCD, are still 
missing to date. 

In Ref. |B] a strategy for evaluating the neutron EDM on the lattice induced by the strong 
CP violating term (|T)) was presented. Such a strategy is based on the standard definition of 
the neutron EDM, involving the insertion of the topological charge {g'^/327r'^) J d^x GG{x) 
in the presence of the charge density operator Jq (see Eq. (4) of Ref. [HI). In case of three 
flavors with non-degenerate masses a complete diagrammatic analysis was performed jO] 
showing how the axial anomaly governs the replacement of the topological charge operator 
with well-defined insertions of the fiavor-singlet pseudoscalar density. The applicability of 
the method to the case of lattice formulations that break explicitly chiral symmetry, like the 
Wilson and Clover actions, was discussed in Ref. using general arguments. 

The aim of this work is twofold: i) to present a complete one-loop calculation of the 
renormalization constants appearing in both singlet and non-singlet axial Ward Identities 
(WPs) using Wilson fermions with the Clover 0(a)-improvement of the action; ii) to 
provide a definition of the gluon operator GG on the lattice outside the chiral limit, taking 
into account its possible power-divergent mixing with the pseudoscalar density. As for the 
one- loop matching, our calculations reproduce all the known results and add: 1) the 

^The present experimental upper limit on the neutron EDM is = \dN\ < 6.3 • 10^^^ (e • cm) at 90% 
confidence level which corresponds to a severe bound on the magnitude of 9. Indeed, using the available 
theoretical estimates from Refs. one has « 3 ■ 10~^^ 1611 (e • cm) leading to \9\ < 2 ■ 10""'^°. The 

smallness of the parameter 9 is usually referred to as the strong CP problem. 
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complete mixing of GG with the divergence of the singlet axial current; and 2) the use of 
both local and extended definitions of the fermionic bilinear operators. As for the calculation 
of the mixing between a lattice discretization of GG and the pseudoscalar density in case 
of lattice formulations breaking chiral symmetry a non-perturbative procedure is presented. 
Finally, as a separate issue, the renormalization properties of the composite insertion of the 
topological charge operator / d'^x GG{x) relevant for the lattice calculation of the neutron 
EDM are discussed. 

The plan of the paper is as follows. In Section |21 the structure of singlet and non-singlet 
axial WPs using the Wilson and Clover lattice QCD formulations is briefly recalled to fix 
notations and basic definitions. In Sections El and 0] the non-singlet and singlet channels are 
considered, respectively. All the matching coefficients with the continuum MS scheme are 
explicitly calculated at one- loop order both with and without the 0(a)-improvement of the 
action. The complete one-loop mixing of the gluon operator GG with the divergence of the 
singlet axial current is evaluated. The issue of the possible power-divergent mixing of GG 
with the pseudoscalar density is addressed and a non-perturbative procedure for evaluating 
the mixing constant is proposed. Moreover, in a separate subsection, the composite insertion 
of the topological charge operator / d^x GG{x) relevant for the lattice calculation of the 
neutron EDM is considered and its renormalization properties are discussed. Section |S] is 
devoted to our conclusions. Finally, all the Feynman rules relevant for our calculations are 
collected in the Appendix. 



2 Axial Ward identities on the lattice 

In the subsequent discussion we will make use of the following definitions and notations. The 
QCD action on the lattice is defined (in Euclidean space) as 

Slqcd = Sp + Sir (2) 

with Su being the pure Yang-Mills component [Zj and Sp the Wilson fermion action 

Sp = a^^{-^^[^/'(a;)(r-7^)f/^(a;)^/'(x + /i) 

+ tP{x + fi){r + -f^)Ul{x)^{x)]+ij{x){mo + ^)ij{x)} , (3) 

where color and fiavor indices are omitted, mo is the (bare) mass matrix, diagonal in fia- 
vor, and the terms proportional to r are necessary to avoid the fermion doubling. The 
improvement of the action is represented by the Clover term 

Sc = '^sw ^-^^i^{x)a^yP^^{x)ip{x) , (4) 

X ii,v 

with P^^ being the usual lattice definition of the field-strength tensor G^^ [S] 

^..W-jiigj^W-c'J) (5) 
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where the sum is over the four plaquettes in the fi-u plane stemming from x and taken in 
the counterclockwise sense (see also Appendix E)) . 

To obtain axial WPs on the lattice, one starts with the usual definition of the vacuum 
expectation value of an operator 0{xi, Xn) 



{0{xi,...,xn)) = ^ I d[G]d[^lj]d[ij] 0{x^,...,xn) e 

Zn J 



(6) 



where the fields are defined only on the nodes of the lattice, Zq is the partition function 
and 5* = Slqcd + Sc- Performing local non-singlet axial rotations over the fermionic fields, 
namely 



ip{x) - 
%l){x) %l){x) 



1 +m"(x)y75 



ip{x) 



1 + ia'*(x)y75 



(7) 



where A" are the usual SU (3) flavor matrices, and taking into account the invariance of the 
measure of integration in Eq. (jH)), one gets 



{O 



6S 



-> = ( ) 

Siia^ix)) 6{ia''{x)) 



with 



((^ Sita'^ix)) ^ = -^imiix)) + Wix){-,moh,^ix)) + (O [X'^(x) + XS(x)]) . (9) 

In Eq. © we have indicated with the backward derivative in the /i-direction with respect 
to X. The non-singlet axial current is given by 



Alix) 



1 



ij{x)Uf,{x)j^—'ilj{x + /i) + h.c. 



(10) 



and the operators X" and are the chiral variations of the Wilson and Clover term in the 
action, respectively. 



X"(x) = -yE[V^(^)^75t/MW(a^ + /^) + V^(a; + /^)^75t/i(a;)V^( 



XI 



A" 



+ (x X - /i) - 4^/'(x) Y75V'(a;) 



XS(x) 



igoar 



csw^i'{x)—l^Ci^yP^u{x)i>{x) . 



(11) 

(12) 
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Similarly, performing flavor- singlet rotations 



one obtains 



{O 



1 + ia;°(x)75 ipi^x 
Tp{x) 1 + ZQ;°(a;)75 



5S 



:) = {- 



SO 



S{ia%x))' '6{iaO{x)) 



) 



(13) 



(14) 



where 



{O 



6S 



6{ia%x)) 



) = -Al{OA^{x))+2{O^Pix) mo 75V^(x)) + {O [xO(x) + (x)! ) . (15) 



The singlet axial current and the operators X° and are obtained from the corre- 
sponding octet operators (fTIUl - lfT^ . respectively, with the naive substitution A"/2 1. 

Equations (jH))-© and (fT^ - (fT3|) are expressed in terms of unrenormalized quantities. 
The operators X*^, X^, X^ and X^ are dimension-5 operators multiplied by one power of 
the lattice spacing; therefore, in any tree-level calculation they vanish identically in the 
continuum limit. 

With the inclusion of quantum corrections the situation changes drastically. In the 
non-singlet channel, the operators X" and X^ mix with the axial current and with the pseu- 
doscalar density: these mixings result in a finite (multiplicative) renormalization constant 
for the current A^, and in an additive renormalization constant for the bare mass tuq, which 
multiplies the pseudoscalar density. The knowledge of such constants allows to identify the 
correct renormalized mass and axial current on the lattice, recalling that Eqs. (jHl)-Q should 
reproduce the corresponding continuum WI in the limit a — > 0. 

In the singlet channel, the mixings of the operator X*^ with the divergence of the (singlet) 
axial current and with the (singlet) pseudoscalar density are in general different from those 
corresponding to X", due to the presence of diagrams involving closed fermion loops (besides 
the ones present also in the non-singlet channel). This could cause in principle the singlet 
renormalized mass to be different from the non-singlet one: if it were so, we would be in 
trouble to identify the continuum limit of the lattice singlet WI, since in the continuum the 
renormalized quark masses in both the octet and the singlet channel are the same. However, 
using general non-perturbative arguments, it has been shown in Ref. Pj that renormalized 
masses are the same also on the lattice, so that, in this respect, a simple correspondence 
between the lattice WPs and the continuum ones can be established. 

The singlet WI on the lattice must reproduce the anomaly, represented by the term 
2NFgQ/32'K'^GG, where GG = 1/2 e^yp„G''^^G'^p^ and is the usual gluon- field strength 
tensor. In a continuum regularization the singlet anomaly is generated by the fact that the 
regulator introduced during the renormalization process is not chiral invariant. Removing 
the regulator leaves a residual contribution in the action proportional to GG. 
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On the lattice, which is a regulator, the action with mo = r = is perfectly chiral 
invariant, but it reproduces no anomaly, since it describes 16 quark species, with opposite 
chiral charges. A possible way to eliminate the 15 spurious fermions in the continuum limit, 
is represented by the Wilson term, which is however a chiral breaking term. The Wilson 
term generates the anomaly. Indeed, in the WI (fT^ - (fTH|) . the chiral variation of the Wilson 
term, namely X^{x), mixes with GG in such a way to correctly reproduce the anomalous 
term present in the continuum [T^. The mixing is independent of r, as long as r 7^ 0. 

Finally, given the anomalous non-conservation of the singlet axial current, the latter 
suffers an infinite (logarithmically divergent) renormalization. Analogous divergences occur 
for the operator GG, so that one needs a mixing between the two operators to end up with 
finite quantities. Furthermore, on the lattice with broken chiral symmetry, one has to take 
into account that the discretization adopted for GG may mix with the pseudoscalar density. 

Since the WPs ((HI)-© and ()14j) - (fT3j) behave differently under renormalization, they will 
be treated separately in the next two Sections. 

3 Non-singlet axial WPs 

In this Section we analyze in detail the one-loop structure of the WI (jHI)-®- For easy of pre- 
sentation we will consider the case of degenerate bare masses mo, which means {A"/2, mo} 
2moA"/2. First we notice that X"' and X^, having the same dimension and the same quantum 
numbers, should have an analogous behavior in perturbation theory. An analysis identical 
to the one made in Refs. fn\ ^ for X" leads to the identification 

X^ix) + XS(a;) = X%x) - 2M(^^)^(x)y75^(x) + (l - zf'^) d.A^ix) , (16) 

where the operator X"(x) collectively refers to dimension-5 operators that vanish in the 
continuum limit for all matrix elements involving elementary operator insertions. Note that 
the renormalization constant Z^^^ refers to the extended definition of the non-singlet axial 
current, given in Eq. (fTUj) . 

The mixing coefficients M^^^) and [1 - Z)^ '] appearing in Eq. ()16|) can be computed 
at the one-loop level simply by considering the correlator {0{x) ip{y)ip{z)) with amputation 
on the external legs, and O given by X" or X^. Its Fourier transform can be pictorially 
represented by the following diagrams (the corresponding Feynman rules are reported in 
Appendix ^ 
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+ 




+ 




(17) 



One gets 



where 



M^NS) 



(1^ 



C 



r (Tq (2r^ , ^ ^ ^ , 4r% ^ ^ 



+ C 



sw 



(A3A5- A4)} 



(19) 



and Eq is the additive mass renormahzation, which appear in the one-loop approximation of 
the quark propagator [see below Eq. (j2Sl)]- 

The one-loop expression for the (finite) multiplicative renormahzation of the axial current 

is 



1 



167r2 



-n (27r)4 



^ -^(5AiA3 + 4Ai-5A2 + A?-13A3 



2A, ' Ai 



+ 4A4) + ^(A3-6Ai + A?) + cw(^;^(4A 



2A1A2 
2 



A, 



+ ^(10AiA3 - 2A1A4 - 2A2A3 - 8A3 - AAl + 6A4 - 4A7) - ^ 

2 2 

^ (9A3A4 - 4A2 - 2A2) + — V(4A3 + - 6A4 + 2A7) 



2A1A2 



+ 



A 



4A, 



(A1A3 - 6A3 + 3A4) + ^(A3 - 2A4))+C5U^— y^(A4 - A3A 

An 



a: 



2r2Ai 



+ 4^^(Ai(-8A3 - 3Al + I2A4 - 4A7 + 6A1A3 - 3A1A4 - A? A3 



- A3A4 + 4A2)]} . 
The symbols Aj stand for the following functions 



(20) 



A3 = \T.^A%) , A4 = Esin^(|)cos^(| 



Ar = Esin^(|)cos^(|) , As = E sin^(|)cos«(|) , (21) 

which for Ai, A7 are identical to the corresponding functions used in Ref. [Hj, while Ag 
is included here for completeness but used only in Section |3] Numerical results for the 
constants C and Z^^'^ at various values of the Wilson parameter r both for csw = and 
for Csw = 1 are reported in Table [T] 



r 


C 


4^^'^ - 1 


C 


Z'i''^ - 1 




csw=0 


csw=0 


csw=l 


csw=l 


0.0 














0.2 


8.09 


-3.62 


5.29 


-2.67 


0.4 


11.41 


-5.82 


6.20 


-4.13 


0.6 


11.25 


-7.00 


5.11 


-4.87 


0.8 


10.43 


-7.89 


3.96 


-5.37 


1.0 


9.64 


-8.66 


3.07 


-5.75 



Table 1: Numerical results for C and [Z)^ — 1] for different values of r and csw- A factor 
g^Cp/lGn^ is understood for [Z^/'^^ - 1] [see Eq. (pUj)]. 



Equations (fT^ -ipU j) . evaluated at csw = 0, coincide with the ones obtained in Ref. jTTj, 
where the Clover term was not included. An explicit check of the correctness of the Clover 
contributions to Eqs. ()19|) and ()20|) can be obtained from the requirement that the WI 
(jHll-Q with O = 4'{y)4'(yZ) should be satisfied. This condition reads as 

A^^{Z^^^A^^{x)^P{y)^P{z)) = 2(mo - M(^^))(^(x)y 75^(0;) ^(y)^(z)) 
- 6{x-y){(^^^,ij{y)jij{z))-6{x-z){^{y)(^^^^^ , (22) 

where we have taken into account Eq. (fTB|) . In momentum space, denoting the amputated 
Green's functions with the insertion of the operators A^y4° and {ipX^^'^^ip /2) by A^^^ and 
A^^, respectively, one gets 

Zr^^lMip'.P) = 2 (mo - M(^^)) A;(p',p) - i,^S-\p) - S-\p'h,^ , (23) 

where p {p') denotes momentum of the incoming (outcoming) quark, and S{p) represents the 
quark propagator. The diagrams corresponding to A^^ are analogous to the ones considered 
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for X"- and X^, reported in Eq. ()17|) . For A;^^ one has to consider only the first diagram, 
being the pseudoscalar current a local operator. After calculating at one loop all the terms 
in Eq. (j23|) . including the Clover contributions, one is able to check Eqs. (fT8|l - (f20|l . The 
explicit one-loop expression for the lattice inverse quark-propagator S~^{p) is 



S-\p) =i^ + mo- -f^Cp(-So + t^J^iip) + moMp)) (24) 



2 

167r 

with So,i,2 given by 



./-, (21)" IA1A2 2Ai 2A2 2A 



+ csw-^(^i^i - A4) + c|^^(A3A5 - A4)} , (25) 
S,(P) = 7.-F„„„,+2|<i..ln(.X(l-.) + „V.(l-.)) 



+ j^(2 - AO + rV.,.(^(A, - 4A3) + 5i-(-9A, + 2A; 
+ 4A3 + 4A5 - 2Aa + 4)) +r^4„.(^^(- A3A4 + 4A|) 

1 2 

+ ^-;^(-4A3A5 + 2A3A6-4A3 + 4A4- A7) 

+ ^(-2Ai A3 + 9A3 - 2A4)) } , (26) 

^2{p) = 4 7£; - Foooo + J dx Inla^ ml{l - x) + p'^x{l - x)^ 
r d'^q r 4Ai - A2 1 , , A2 

;^)+csw^(-A3A5 + A4) + 4^^^-^(A3A5 - A4) 

2 2 12 

-(-A3A5 + A4)} (27) 



+ 


16n^ 


+ 


r\- 


+ 





Ai 

with Foooo = 4.36898, Foooi = 1.31096 and 7^ = 0.577216. Equations (12323) coincide with 
the analytical results of Refs. [TT1|T2] (at csw = 0) and Ref. ^ (at csw = 1)- 

The check of Eqs. (|18 p -(|20 p via the WI (j23p can proceed now through the projection of 
all the quantities onto the structures {]^' — y)75 and 75 (orthogonal to each other) with the 
appropriate normalizations. In other words we calculate the quantities 

nr(p',p) = ^rr(T(p',p)r) , (28) 
9 



where the trace is taken over the Dirac indices, T{p',p) stands for any of the terms appearing 
in Eq. ()23|). F is equal to — ]i)75 or 75 and A^r is a normahzation constant. We have 
calculated all the relevant projections and checked positively the correctness of Eqs. (|T9|) - 

Note that the projection of the amputated diagram of the (extended) axial current over 
i"^' — and of the pseudoscalar density over 75 can be compared with the corresponding 
quantities calculated in the continuum M5'-scheme in order to extract the matching constants 
between MS and the lattice [see Subsect. ()3.1|) ]. because the dependence on the external 
momenta is cancelled out in the difference between the two schemes. 

Finally, from Eq. (j22j) it can be noticed that the explicit contribution of the operators 
X" and has been traded with a redefinition of the lattice axial current and of the bare 
mass mo, i.e. 

mo — ^ mo - M^^^) = m^ . (29) 



3.1 Matching with the continuum in the non-singlet channel 

The Green's functions calculated with the (bare) fields and operators on the lattice can be 
matched to the corresponding ones renormalized in the MS scheme, by rescaling the fields 

and the operators on the lattice with appropriate constants. For instance the bare quark 

1 

field on the lattice ^ipLix) can be rescaled through ipL^x) = Z^ipR^x), where the subscript R 
indicates a renormalized (MS) quantity in the continuum. The constant can be evaluated 
via the relation 



{s-\p))^ = z4S'\p)Y 



(30) 



which simply follows from the definition of the quark propagator. At one-loop order the 
relations 



(5-^(p))^ = ¥ 
{S-\p))^ = t^ 



16n' 



imply 



1 - 



_9o_ 

IGtt' 



+7711 



-rriR 



As, = Sf (p) - S[(p) . 



1 - 



167r2 



CpEf(p) 



(31) 



(32) 



Analogously, one can define a matching for the mass on the lattice from the relation 

rriR = ZmrriL, (33) 
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obtaining, thanks to the use of Eq. ()3H) and of the definition 

Z„ = Z^l - JI^C^AsJ . (34) 

A similar procedure leads to the definition of the Z's for the matching of composite op- 
erators. In the non-singlet channel, we are interested in Green's functions involving fermion 
bilinear operators O contracted with two external quark fields, with amputation on the ex- 
ternal legs. The corresponding functions in momentum space will be indicated by Ao{p',p)- 
One gets the relation 

A^{p',p) = Z^ZoAUp',p) (35) 

which defines Zo as the constant of matching for the operator O. In perturbation theory 
one can express the constants in Eq. as 

Z^Za = 1 + t^CfAo , Ao = Ag - , (36) 
where the expressions A collect all the one-loop diagrams. The relation (jHUj) implies 



Zo = 1 + t^Cf{Ao + A^,) . (37) 



167r2 

Denoting the subtraction scale in the continuum by the symbol /i, our results for A^^ 
and Asj are 

d^q ( As I A3 



A,, = l-ln«V-7i. + i^oooi-16vr^£7^{^ + 



(27r)4ll6A? 8A1 8A1A 



3 

2 



+ ^^2:^(2^4 - A3A5) + -—(2 - Ai) + r'cswiTTin^i^^ " ^As) 



SAfAa' 4A2' ' "''^V2AiA 

+ -^(-9Ai + 2A2 + 4A3 + 4A5-2A6 + 4) 

+ ^'^^-(sAk^"^^^^ ^ '^'^ + 8Ak^"'^^^^ + 

- 4A3 + 4A4- A7) + ^(-2AiA3 + 9A3-2A4))}, (38) 



A 



S2 



/ 1 , 2 2 \ 2 r d'^Q ( 4Ai - A2 1 



(27r)4 I 4Af A2 4A2 
+ ^(2A,A. + ^-4A3)+r*(-2|) 

+ cw-r2(-A3A5 + A4) + c|vK , . . (A3A5-A4) + 
A2 4A1A2 

+ 4^y^p(-A3A5 + A4)} , (39) 
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The numerical results obtained for Eqs. (jHHj) and (jSH) at various values of r, are 

reported in Tables El and El for the two cases csw = and csw = 1; respectively. 



r 


So 






0.0 





-6.04 


33.17 


0.2 


-19.79 


-7.13 


18.63 


0.4 


-30.70 


-8.91 


8.35 


0.6 


-38.29 


-10.47 


3.86 


0.8 


-44.96 


-11.77 


1.55 


1.0 


-51.43 


-12.85 


0.10 



Table 2: Numerical results for Sq and A^i j [see Eqs. (j^ . (p?^ and (jH^jl ] for different values 
of r and with csiy = 0. The choice = 1/a is understood. 



r 


So 




Ae. 


0.0 





-6.04 


33.17 


0.2 


-12.04 


-6.52 


21.81 


0.4 


-18.31 


-7.40 


14.69 


0.6 


-22.98 


-8.16 


11.95 


0.8 


-27.44 


-8.75 


10.74 


1.0 


-31.99 


-9.21 


10.10 



Table 3: The same as in Table |21 but with csw = 1- 



As for the qi 
A, = 

+ 

+ 



+ 



lantities A/, A^^ and A^^^^^^ one obtains 



^7: 

4{-l + In - i^oooi + 1e - levr^ J 



2 r A 



A, 



A.A 



3^5 



(27r)4 L i6Af 4AiAi 



A? 



4A 



A1A5 + 3A3) + r^—^ + (A3A5 - A4 



r^A. 



4A 

A3A5 + A4) + c 



A 



2 

sw 



r^Ai 

4Ar 



"^^iA^Af 
4{-l + InaV^ - i^oooi + 7b - IStt^ 



(A3A5-A4)]} 



A3 ^ A3A5 



-n (27r)4 L 16A3 4AiA2 



4A; 



r(AiA5 + A3)+r 



.A^ 
4Ai 



2 ^^Ag 

+ ''^^4A;aI 



(A3A5 - A4) 



r^Ai 
4A^ 



(A3A5-A4)]} 



(40) 



(41) 
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Kext 



-2 + In - i^oooi + 7b - IStt^ /" 

J— 7r 



ci^g r A3 



A, 



(27r)4L l6Af AiAi 8A1 



+ 2A4)+r^-^ + 



Ai 

4A^^ 2A|' 



-3Ai A3 - 4Ai A5 + 2A1A6 + 4Ai + 3Al 



^(2 - Ai) + C5H/^(-7AiA3 + 4A1A4 



Al + 7A3 - 2A4) 



+ 2A2A3 + 8A3A5 - 4A3A6 - 4A3 - I5A4 + 6Av 



+ ''^^2AiAi 



(3A3A4 - 4A3 Ay - 4A^A5 + 2A^A6 + iAj + 2A\ 



+ Ai(3A2-A3A4- AiA2)) + c; 



sw 



r^Ai 

2A^ 



(4A3A5 - 2A3Ae 



I2A4 + 4A7 - 4Ai A3 + 3Ai A4 + A? A3 



y(NS) 



We notice that using Eqs. (j5H|) and (j42|) one has 



167r2 



7m 75 ' ^1 



- 1 



(42) 



(43) 



where Z^^"^'* is given by Eq. (j^ . 

For completeness, we report also the matching constants for the operators ■?/'(x)7^'?/'(x), 
^(3^)7^75^(3;) and ^{x)af,^'4){x) [H]: 



A^^ = -2 + In a V' - i^oooi +1e- leTr^ T 

J — IT 



- d\ r A3 



A, 



(27r)4L I6A? AiA^ 



+ 
+ 



4^? , r'^csw 



2AI 



(A1A5 + 3A3) + + (-A3A5 + A4) 



(A3A5 - A4) 



'7m 75 



2AiAi 

-2 + InaV^ - ^0001 + 1E- 16vr2 T 



(44) 



- rf^g r A3 ^ A4 



(27r)4L i6Af AiA^ 



+ 
+ 



I (-A1A5 + 5A3) + r^^ + ^ (A3A5 - A4) 



2AiAi 



+ 



-167r 



- (A4 - A3A5) 
d'^q V 1 



(45) 



AiA3-4(A3-A4) +2r^-| + r^ 

^2 ^2 



(27r)4 L3AiAi 

(A4 - A3A5) + ''!ffr^ (A4 - A3A5 



2 A3 



.Al 



^^4h/A i 
3Ai ^""^ • 3Ai 



(46) 
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The numerical values of the matching constants A/, A^^, A^^, A^^^g, A^r^^ and A!^^^^^, 
obtained for various values of r, are collected in Tables |^ and El Our results for the local 
operators have been explicitly checked against the corresponding ones of Refs. [HI El EI, 
while those for A'^^t^^ at csw = coincide with the ones of Ref . . 



r 


A, 


A,, 








J\ext 

7f«75 


0.0 


-33.16 


-33.16 


-8.74 


-8.74 


0.74 


6.04 


0.2 


-18.63 


-26.72 


-8.97 


-4.92 


-0.37 


3.52 


0.4 


-8.34 


-19.75 


-8.76 


-3.05 


-1.86 


3.09 


0.6 


-3.85 


-15.10 


-8.36 


-2.73 


-2.90 


3.47 


0.8 


-1.55 


-11.98 


-8.02 


-2.80 


-3.63 


3.88 


1.0 


-0.10 


-9.74 


-7.76 


-2.94 


-4.17 


4.19 



Table 4: Numerical results for various A^ for different values of r and with csw = 0. 



r 


A, 








A, 


/\ext 

7u75 


0.0 


-33.16 


-33.16 


-8.74 


-8.74 


0.74 


6.04 


0.2 


-21.80 


-27.09 


-8.27 


-5.62 


0.23 


3.85 


0.4 


-14.69 


-20.89 


-7.45 


-4.35 


-0.61 


3.27 


0.6 


-11.94 


-17.05 


-6.82 


-4.27 


-1.23 


3.29 


0.8 


-10.73 


-14.70 


-6.40 


-4.42 


-1.64 


3.38 


1.0 


-10.10 


-13.17 


-6.12 


-4.59 


-1.93 


3.46 



Table 5: The same as in Table EJbut with csw = 1- 



In order to get a further independent check of the expression ^T^i for the constant C one 
can consider the relation 

mR = Zmimo- -^Cf — = , (47) 

y 167r^ a J z^^ 

where the first equality follows from the definition (jH^ . while the second is a consequence 
of Eq. ()22j) and of the definition (jSHJ. Now, using the one- loop expressions 

Zm = 1--|^C^(As, + As/ 
167r^ a 



^75 = 1 + t|^C7^(A,,+As,) , (48) 



1677 
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one easily finds 

C = -AS2 - A75 , (49) 
which can be checked analytically using Eqs. (fT^. (jH^ and (jUJ. 

4 Singlet axial VF/'s 

Let us now turn to the singlet WI (fT^ - (fT3j) . Because of the anomaly, the operator X° + Xfl, 
representing the chiral variation of the Wilson and Clover terms, can be written as ^T] 

X°(a;)+X°(x) = X\x)+X''c{x)-2M^^^{x)^,,p{x) + {l-Zf)d^A^{x) 

+ 2Nf^Z^^GGs^,{x) , (50) 

where, again, X° and indicate dimension-5 operators vanishing in the continuum limit in 
all matrix elements involving elementary operator insertions, and the symbol GGsub indicates 

a lattice discretization of the corresponding continuum operator GG. In Eq. M is the 

(s) 

additive mass-renormalization constant, while Z\ and Z^^ are multiplicative constants in 
the singlet channel. 

Before addressing the calculation of the renormalization constants, we point out that in 
Eq. (I^Uj) GGsub should be defined in such a way to avoid ambiguities with the pseudoscalar 
density. Indeed, let us consider the commonly used discretization GG^ obtained through 
the use of the symmetric plaquette, viz. [TB] 

GGiix) =e^yp„Tr{P^^{x)Pp„{x)) , (51) 

where P^i, is defined in Eq. © and the sum over the greek indices is understood. Expanding 
the definition in powers of a, one gets 

GGl{x) = GG{x) + a^Oeix) (52) 

where GG = l/2£^j,po-G'^j,G^o. and Oq denotes a dimension-6 operator. Since with Wilson 
fermions chiral symmetry is broken, there is no guarantee that the operator Oq{x) does not 
mix with the pseudoscalar density. In case of mixing one expects by dimensional reasons 
a (cubic) divergence, leading to a linear-divergent mixing of GGl with the pseudoscalar 
density. Such a divergence has to be subtracted and therefore the operator GGsub, appearing 
in Eq. (jHUjl . should be related to GGl by 

GGsubix) = GGLix) - AMz.V^(x)75^(x) (53) 

where AMl contains the 1/a divergence and is in general dependent on the particular 
discretization GGl- Up to now, a non-perturbative definition of AMl is not known, and 
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we anticipate here that at one-loop order we find no mixing, i.e. AM^ = OIqq). However a 
non-perturbative definition of AMl is mandatory, particularly outside the chiral limit. To 
this end we require that GGgub should be a chiral singlet. 

First of all our requirement on GGsub is needed to guarantee the proportionality between 
[mo - M^^")) and (mo - M^*^^), as it is demonstrated in Ref. p. This implies that: i) the 
chiral point is the same in the singlet and non-singlet channels; and ii) the renormalized 
mass is the same in the two channels and it can be therefore identified with the continuum 
renormalized mass. The latter point follows from a suitable choice of the renormalization 
constants for the singlet and octet pseudoscalar densities, which in turn arises from the 
transformation properties of both scalar and pseudoscalar densities as members of the same 
chiral multiplet. 

Let us now consider the non-singlet WI and the case O = S"'{yi,y2) GGsub{z), where 
S°'{yi,y2) = V^(?/i)A"V^(y2)/2 with yi ^ yi\ one has 



Al{^Al{x)S\y^,y2) GG.^,{z)) 



= 2mL{P^{x)S\yi,y2)GGsub{z)) 

- Cf [5{x - yi) + 5{x - y2)\ {P^iyi, I/2) GGsub{z)) 

+ +X^(x)l 5^(1/1,1/2) GGsubiz)), 



(54) 



where P^{x) 

K7VS) 



^P{x)\^-f5^{x)/2, P5(z/i,Z/2) = ^(1/1)75^(1/2), Alix) 



rriL = 

mo — M^""^ and the sum over the fiavor index a (a = 1, ...,8) is understood. In Eq. (j^lj) 
we have considered that 6GG subiz) / S^ia"" (x)) = 0, because GGsub is assumed to be a chiral 
singlet. After integration over the whole space-time in x and one gets 



= m^ij d'x d'z P^{x) 5^(1/1,1/2) GGsub{z)) - CF{P,{y,,y2) J d^z GG^ubiz)) 



+ ^{jd^xd'z X''{x)+X^{x) S^{y,,y2) GGsubiz)) . 



(55) 



We have now to discuss the possible presence of contact terms which may arise in the first and 
third terms of the r.h.s. of Eq. when x ^ z. Note that the operator S°'{yi,y2) is a string 
of elementary operators taken at different space-time points yi 7^ 1/2; therefore it cannot 
generate contact terms when inserted with composite operators. The general structure of 
the possible contact terms in Eq. can be derived using the results of Ref. where 
the method of functional integral with the generalized mass term, a°''il)\°"ip /2 + m°"il)\°''^^il) /2, 
was developed to derive VTJ's. The only possible contact terms should be as follows 



j d^x dhP^{x) ^'^(1/1,1/2) 0'{z) 



d X d z 



X\x)+X:^{x)\ S\y,,y2) 0\z 



'^The integration over x is used to cancel out the l.h.s. of Eq. H54() . including possible contact terms 
between v4^(a;) and GGs„b(z). The integration over z allows to get rid of the mixing of GGs„b(z) with 
the singlet ^^iA^^(x) [see below Eq. (|59|l ]. Note that such a mixing does not change the property of the 
renormalized GG^ to be a chiral singlet. 
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J d'x S^{yuy2) St Ax) , (56) 



where 6t7c^{x) is the chiral variation of the mass field 7r''(x) and 0'{z) is an operator with 
the same transformation properties and dimension of GGsubi^), i.e. a pseudo-scalar, chiral- 
singlet dimension-4 operator, build up with the fields ip{z), ip{z), a°'{z) and 7r"(z) [T7]. 
Since 0'{z) should be a chiral singlet, the third form of the contact terms (j56p is vanishing. 
Moreover, the only dimension-4, chiral singlet operator allowed is just the generalized mass 
term, i.e.: O' = a°'iljX°"ip/2 + z7r"'0'^"75'0/2, which however is not pseudoscalar. Thus, we 
conclude that no contact terms are present in Eq. ()55p. Finally, according to Refs. fn\ ITT] 
in the continuum limit the third term in the l.h.s. of Eq. ()55|) vanishes and therefore it will 
be disregarded. 

Substituting GGsub{z), defined in Eq. into Eq. one gets 

-GFiJd'z P,{yi,y2) GGL{z))+mL{Id'x d'z P^{x) ^"(^1,^/2) GGl{z)) 
^ -CFUd^zP,{y,,y,)P,{z))+mLUd^xd^zPi{x)S-{y,,y2)P,{z)) ' ^ ^ 

Note that our non-perturbative definition of AM^ is gauge- invariant because the loss of gauge 
invariance is given by a factor which is the same in the numerator and in the denominator 
of the r.h.s. of Eq. 

Let us now address the calculation of the renormalization constants M , Zj^ and Z^g 
relevant to the singlet WI. Substituting Eq. (jHOJ into the WI (fTij) - (fT^ . one gets 

A^(zfA^(x)0) = 2(mo-M(^))(V^(x)75V^(a:)C») 

+ 2iV.ji,(Z,gGa„.(x)0)-(^^), (5^ 



where we have dropped the term involving X° + Xq. Because of its anomalous non- 
conservation the singlet current Afj_ suffers a logarithmically divergent renormalization [18\ . 
A mixing between Z^^A^^ and Z^g GGgub is needed to obtain finite operators jU] 



A^{x) = zi'^A,{x)~{glZc)zi'^A,{x 



2Nf^GGr{x) = 2Np^Z^QGGsub{x)-ZcZ)^>d^A^{x) , (59) 

where we stress that GGsub{x) = GGl{x) — /S.M Lip{x)^^ijj{x). The above redefinition can be 
performed into Eq. (j^Uj) by adding and subtracting the counterterm ZcZ^^'^ ■ d^A^{x) 0. 

We have evaluated at one loop level all the constants appearing in Eq. (jSHI). The con- 
stant Z^^-* can be obtained by evaluating (X°(x) + X^{x)il){y)il){z)): to 0{gl) its value is 
identical to that of Z^^\ Eq. ()20|) . while beyond the one-loop order the two constants will 
differ because of the contributions (in the singlet channel) coming from fermion loops. The 
same is true for M^^\ which therefore at one- loop is given by Eq. ()18|1 . as well as for the 
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renormalization constants of the pseudoscalar density operators in both the singlet and non- 
singlet channels. This implies that at one- loop the renormalized masses (mo — M )/Z'^f' 

and (mo — M*'^'^^)/Z^^'^^ coincide and match the corresponding renormalized mass in the 
continuum. We stress that the equality between the renormalized masses in the singlet and 
non-singlet channels is a more general result, valid at any order of the perturbation theory 
and resulting from the property of GGsub to be a chiral singlet, as illustrated in Ref. 

The mixing between X°(x) + Xq{x) and GGsub{x) can be computed at one loop by 
evaluating the correlator ([X°(x) + X^{x)] G^p{y)G^:^{z)) . This was carried out in Ref. [lOj . 
where it was shown that the apparent dependence of the correlator on r actually disappears 
as far as r 7^ 0, and that the tree level expression 2NpgQ/32'7i'^ GG is reproduced exactly. 
Thus, from Eq. (jKHjl one has Z^g = 1 + O^Qq), and, since the logarithmic divergence in the 
singlet axial current manifests itself at the two-loop level, it follows that QoZq = 0{gQ) and 
consequently Zq = 0{gl). 

The same considerations hold as well for the matching constants between the continuum 
and the lattice. For Z^g and Zc appearing in Eq. we can write 

Zgg = l+^?o4?+0(^?o') , 
Zc = J^C^(45ln(aV) + 45)+0(^?o'). (60) 

The coefficient z^^J can be computed at one-loop in a pure gauge theory (see Ref. [12] and 
references therein quoted), and at this order it is unaffected by the presence of the Clover 
term in the action"^. One obtains jTH] 

4? = ^(-4^ + ^0000 + ^ + ^) (61) 

where N is the dimension of the gauge group and Zqooo = 0.15493. The constants z^^ and 

~(2) 

Zg^ can be computed by evaluating, in the continuum and on the lattice, the correlator 
{GG{x)ilj{y)ip{z)) , with amputated external quark propagators. On the lattice GG should 
be GGsub given by Eq. (j^ . However, we have explicitly checked that the operator GGl 
does not mix at one-loop with the pseudoscalar density and, therefore, in our leading-order 
approximation there is no difference between GGsub and GG^. Thus, on the lattice we have 
to compute the following diagrams: 



FT{GGL{x)iP{y)^{z)) 




■^Beyond 0((7q) the presence of the Clover term in the action will, in general, affect also Z^~. 
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where the ijji/jg and ip'ipgg vertices include the contributions of both the Wilson and Clover 
action. The second diagram gives a null contribution, since the insertion of GG contracted 
with two gluon lines is antisymmetric in the Lorentz indices of the gluons (via the e-tensor) 
and symmetric in their color indices, whereas the ipipgg is symmetric in the Lorentz indices 
for the Wilson case and antisymmetric in the color matrices for the Clover case (see the 
Feynman rules in the Appendix). The first diagram displays a log divergence as a —> 0, in 
analogy with the continuum calculation. The difference between the two schemes gives the 
following result 



,(2) 



- 6 



~(2) 



+ 
+ 



+ 



22 -6{^E- Foooi] 
1 



1 

4A1A2 



(1OA3A4 



+ r csw 



(45A3 
1 



6A^ 



- IGtt 
3OA3 
4OA4 



f d^q 


r3A3 


f 1 M 


1 (27r)4 


^A? 


UAi A2/ 



19A^ + 64A4 - 46A7 + 12Aj 



12A- 



4A, 



[2A1A3-I7A3 + 6A4) 



2A1A2 
2A4A6 - 8A5A7 
1 



8A3A4 - 3A3A^ + 3A3A6 - 4A^A5 + IOA4A5 + 2A4A^ 



I4A7 + I2A8 + 



A, 



I 22 



4AfA2 



IOA3A4A 



2A3A4A^ 



(A3A2 - A3A6 - 2A4A5 + 2A, 
2A3A4A6 + 8A3A5A7 



+ I4A3A7 - I2A3A8 - 8A^A4 + BA^A^ - SA^Ae + 4A^A5) 
^ (2A3A4A5 - 2A3A7 - A^A^ + A^Ae' 



2A1A2 ^ 



(2) 

The numerical values of the constant obtained for various values of r at csw = and 
Csw = 1 are reported in Table IHl 



(63) 



4.1 Composite insertions of the topological charge and the neu- 
tron EDM 

In this subsection we take the opportunity to briefly address a separate issue concerning the 
composite insertion of the topological charge relevant for a lattice evaluation of the neutron 
EDM induced by the strong ^-term (^. Following Ref. [Hj the standard definition of the 
neutron EDM, d^, involves the insertion of the topological charge ((7^/3271^) / d^x GG{x) 
in the presence of the charge density operator Jq. Treating the ^-term as a perturbation at 
first order, one has 



dpj = —iO- 



3_ 

327r2 



d'yyo{N\Uy) 



d^x GG{x) 



\N), 



(64) 
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r 




^(2) 




csw=0 




0.0 


19.82 


19.82 


0.2 


19.72 


19.93 


0.4 


19.11 


19.98 


0.6 


18.11 


19.94 


0.8 


16.95 


19.87 


1.0 


15.77 


19.81 



Table 6: Numerical results for z^J, [see Eq. (jHSI)] for different values of r and csw- 

where |A^)o is a shorthand for \N)g=Q. In case of three flavors with non-degenerate masses 
a complete diagrammatic analysis was performed in Ref. jO] showing how the axial anomaly 
governs the replacement of the topological charge operator with well-defined insertions of 
the flavor- singlet pseudoscalar density*^. 

Thus the question is the possible presence of contact terms in the composite insertion 
of the topological charge operator with the electromagnetic (e.m.) current operator, which 
would lead to ambiguities in the numerical evaluation of Eq. ()64j) . The contact terms in 
Eq. ()64|1 should be of the form OI^^q{x) ■ 5^^^ (x — y), where is a local operator of dimension- 
3, which transforms as a non-singlet pseudo-vector and is conserved. The last property 
derives from the fact that the e.m. current is conserved both with and without the ^-term 
in the QCD action at any value of the parameter 6. The non-singlet nature of Oj^ is related 
to the non-singlet nature of the e.m. current operator. The only candidate for Oj^ is the 
non-singlet axial current which however is not conserved (outside the chiral limit). 

The absence of contact terms in Eq. ()64|1 can be derived also using the generalized mass 
insertion method of Ref. J7] applied to a vector WI with the operator O given by GG. 

5 Conclusions 

A complete one-loop calculation of the renormalization constants appearing in both sin- 
glet and non-singlet axial Ward identities using Wilson fermions with the Clover 0(a)- 
improvement of the action has been performed. Our calculations include: 1) the contri- 
butions arising from the Clover term of the action; 2) the complete mixing of the gluon 
operator GG with the divergence of the singlet axial current; 3) the use of both local and 
extended definitions of the fermionic bilinear operators. 

In the singlet channel a definition of the gluon operator GG on the lattice outside the 
chiral limit has been proposed. Our definition takes into account the possible power-divergent 

'^We point out that the final result of Ref. [S] crucially depends on the equality between the renormalized 
masses in the singlet and non-singlet WPs. We stress once more that such an equality follows from the 
property of GGsub [see Eq. l|5^ ] to be a chiral singlet (cf. Ref. 
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mixing with the pseudoscalar density, generated by the breaking of chiral symmetry. No 
mixing has been found at one-loop order and a non-perturbative definition of the mixing 
constant has been developed. 

Finally, the renormalization properties of the composite insertion of the topological charge 
operator / d'^x GG{x) relevant for the lattice calculation of the neutron electric dipole mo- 
ment have been discussed, showing that no contact terms arise when the topological charge 
is inserted with a conserved current. 
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Appendix: Feynman rules 

The Wilson action is given by 

+ + /i)(r + 7^)f/^(a;)^/'(x)] + ^(a;)(mo + , (65) 

with the following notation {7^,7,.} = 25^u, a^u = Yl^i.^u]/'^, U^{x) = exp[igoaGf,{x)], 
G^ = Gp'' and Tr(rt'') = 5ab/2. 

The action is improved to 0(a) via the Clover term [201 

Sc = -a^^Yl '^sw '^-^^i^{x)a^,uP^y{x)ip{x) , (66) 

X 11,11 ^ 

with P^y being the usual lattice definition of the field-strength tensor G^y ^ 

where the sum is over the four plaquettes in the /i-z/ plane stemming from x and taken in 
the counterclockwise sense, i.e. 

Ui = U^x) Uy{x + 12) Ulix + u) Ul{x) , 

U2 = Uy{x) Ul{x - fi + u) Ul{x - jj) U^{x - /i) , 

Us = Ul{x - fi) Ul{x - fi - u) Uf,{x - fi - u) U,,{x - u) , 

t/4 = ulix - v) U^[x - v) U,[x + /X - I/) Ul[x) . (68) 
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Propagators 



5^^ -(!-«/ ^2 



with 



a 2 



"2 \ " -2 



• »> • 



[~ sin(pAa) + ("^o + ~" ^^^^ ' 



QCD vertices 

The indices W and C refer to the Wilson and Clover action respectively 




-go 



{i7/. cos((g + q')f,^) +r sin((g + } 
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2 

= ^^palt", t%j [n, sin((p + ) -r cos((p + p')^) 



Extended operators 

For simplicity we report only the Feynman rules for flavor-singlet operators. 
Operator A^Ajj^^x): 



aAf,Af,{x) = - Af,{x - /i) 

1 



^m(^) = + i^hi^i5Ul{x)ij{x) + 'ii){x)^^^^Uf,{x)'ii){x + ^J) 



= -- 1^7/^75 (sin(p^a) -sin(p^a)) , 
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Operator X^: 



2a 

+ {x ^ X — — 4:ll!{x)j5lp{x) 



■ ► 



-75 - cos(p^a) - cos{p'i^a)) , 




igor 



—ippa 1 
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g^ar 



-75 



[t"^ t%j H { [^pa(^pu sin((g - g')pf ) (sin(gi.a) + sin(g») 



+ 2^P. 



cos((gp - g^)-)-cos((gp - g^)- - q'^a 

- cos((gp - g;)| + g^a)-cos((gp - g^)^ + {q„ - q'p)a 
+ cos((gp + g;)^)-cos((gp + g;)| -g;,a 

- cos((gp + g;)|-g^a)-cos((gp + g;)|-(g<, + g|,)a)]} . 



(86) 



Operator GGl- 



GGl — SfjLupaTr{PpyP^ 



per) 



(87) 



2^ab^£uo,fii3COs{py-)cos{rp,-)si\i{paa)sm.{rpa) . (88) 

"a/3 II 
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